notebook.community

Edit and run

Running a sanity test of my mock testing pipeline by putting an HOD on the MDPL2 pipeline and seeing if i get the same thing out as my emulator predicts.



In [1]:

    
import numpy as np
from pearce.mocks.kittens import MDPL2









    



/home/users/swmclau2/.local/lib/python2.7/site-packages/h5py/__init__.py:36: FutureWarning: Conversion of the second argument of issubdtype from `float` to `np.floating` is deprecated. In future, it will be treated as `np.float64 == np.dtype(float).type`.
  from ._conv import register_converters as _register_converters



In [2]:

    
cat = MDPL2()



In [3]:

    
cat.load(1.0, HOD='zheng07')



In [4]:

    
true_hod_params = {'alpha': 1.083, 'logM0': 13.2, 'logM1': 14.2, 'sigma_logM': 0.2}



In [5]:

    
true_hod_params['logMmin'] = 13.0 # TODO fit nd
nd = 5e-4

cat.populate(true_hod_params)

print len(cat.halocat.halo_table)

welp, gotta move this to a script, keep killing the kernel trying to load things up.